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^ 1 Overview 

. The talk given by the author at the CRM workshop on Superintegrability in Sep. 2002 and this related paper 
|Jh . report on work in two subjects. One is the collaboration with Vladimir Sokolov and Takayuki Tsuchida in 
i ^, 1 an effort to classify polynomial integrable vector evolution equations. The other is the computer algebra 
package Crack which did the main computations in solving large bi-linear algebraic systems. Although 
originally designed to solve over-determined systems of partial differential equations a number of extensions 
CN ■ made Crack a powerful tool for solving systems of bi-linear algebraic equations. Such systems turn up in 
many different classification problems some of which were investigated by other participants of this workshop. 
In sections El and El two additional applications are outlined. 

In the talk on which this article is based a method to reduce the length of equations was presented which 
proved to be useful in solving the bi-linear algebraic systems. Due to numerous asked questions about the 
computer program, a more complete overview is given in the appendix. 
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' 2 The classification of integrable vector evolution equations 

>< 

?h ' The method to use symmetries to classify non-linear evolutionary 1+1 dimensional PDEs is the most pro- 
ductive one and known for some time (see [U El El E] ) • 

An extension of the simplest approach E] to the case of so-called vector evolution equations is de- 
scribed in work with Vladimir Sokolov jJJJ. Examples of vector evolution equations are two different vector 
generalisations 

U t = U XXX + (U,U)U X , (1) 
U t = U XXX + (U,U)U X + (U,U X )U (2) 

of the mKdV equation where U(t,x) is an iV-component vector and (•,■) stands for the standard scalar 
product. 

In the performed symmetry classification we considered equations of the form 

U t = f n U n + /„_! U n -! + ■ ■ ■ + h Ux + / U, Ut = (3) 



1 



where U = (U 1 , U 2 , . . . , U ) is an unknown vector of arbitrary dimension N and coefficients fi are polyno- 
mials of scalar products (Ui, Uj) , < i < j < n. 

For more details on possible scalar products, orthogonal symmetry transformations and generality of iV 
see jnj. 

All such equations were determined that omitted a symmetry of the same form 

U T = g m U m + g m _ x U m -i H \-9iUx + g U 

where the differential order n of the equation and the order m of the symmetry have selected values. By 
taking fa, gj to be homogeneous polynomials of the scalar products one achieves that the symmetry condition 
Ut T = U T t yields an over-determined system of bi-linear algebraic conditions for the un-determined coefficients 
of both polynomials f\ and gi. 

The number of coefficients of the equation and symmetry is reduced further by a homogeneity assumption 
that the differential equations and symmetries are invariant under the scaling group 

(x, t, U) — ► (a _1 x, a~ M t, a x U). 

with one value A for all components of the vector U. 

Finally, we restrict A based on results for the scalar case proven in [H] that a A-homogeneous polynomial 
equation with A > may possess a homogeneous polynomial higher symmetry only if A G {2, 1, |}. 

The differential orders of the equation and symmetry that have been investigated are also motivated by 
the scalar case where infinite commutative hierarchies have either a lowest order of the equation of 2 and 
symmetry of order 3 (in the following called type (2,3) like the Burgers equation) or an equation of order 3 
with a symmetry of order 5 (type (3,5) like the Korteweg-de Vries (KdV) equation) or an equation of order 
5 with a symmetry of order 7 (type (5,7) like the Kaup-Kupershmidt and Sawada-Kotera equations). 

For the nine resulting cases shown in table 1 computer programs generated the equations and symmetries, 
computed the commutator and formulated and solved the bi-linear system automatically (apart from the 
largest case A = |, type (5,7) where the solution was not fully automatic). 
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Table 1. A classification of single vector equations. 



Comments: Times are measured on a 1.7GHz Pentium 4 running a 120MByte Reduce session under 
Linux. Nonlinear A = 2 equations of order 2 do not exist. All found solutions have already been known: 
equations (JTJ), (J2J) are vector generalizations of the mKdV equation, and equation 

U t = U xxx + 3(U, U)U XX + 6<Z7, U X )U X + 3(U, U) 2 U X + 3(U X , U X )U (4) 

is a vector analogue (reported in [TU]) of the Ibragimov-Shabat-Calogero equation [ZHHj- Solutions of order 
(5,7) are only symmetries of found lower order equations and have therefore not been listed in the table. 
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3 NLS Systems with two Vector Unknowns 

More successful has been the investigation of integrable vector NLS-type systems of the form 

U t = Uxx + PiU x + p 2 V x + p 3 U + p 4 V, 
v t = -v xx + p 5 U x + p G V x + p 7 U + p 8 V, 

where U and V are vectors and the coefficients Pi are A-homogeneous polynomials depending on all possible 
scalar products of vectors U,V,U X ,V X . For A = 2 second order systems can only be linear but for A = 1 or 
1/2 results are shown in table 2. Just as in the scalar case (see [2]) a symmetry of the form 

U T = U xxx + qiU xx + q 2 V xx + q 3 U x + q 4 V x + q 5 U + q 6 V, 

(6) 

V T = V xxx + q 7 U xx + q 8 V xx + q 9 U x + q 10 V x + q u U + q 12 V, 

is assumed where the coefficients qi are A-homogeneous polynomials of all possible scalar products of 

U) V : U X i V x ^ U XX i V xx - 



A 


1 


1/2 


orders of (cq,sym) 


(2,3) 


(2,3) 


# of unknowns (eq,sym,tot) 


9,15,24 


53,155,208 


# of equations 
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total # of terms 


242 


28768 


av. # of terms/equ. 
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23.8 


time to formulate 


2.3s 


lh 


time to solve 


2.9s 


26m 15s 


# of solutions 
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Table 2. A classification of NLS-type systems of vector equations. 



Results: 

The 2 solutions for A = 1 are known (jll^H2|). They are special cases (see [3]) of a generalization of the NLS 
system by Svinolupov using a Jordan triple system (in [T3*j). For A = | after identifying solutions through 
U <-> V, t <-> —t six solutions remain. Two of these systems had been known (|22j). For the remaining four 
systems, Takayuki Tsuchida showed C-integrability for two of them and S-integrability for the other two. 



4 Systems with one Scalar and one Vector Unknown 

Computations to classify single vector equations that involved an arbitrary constant vector did not give 
new results but the possibility to apply an orthogonal transformation to make the constant vector equal 
(1,0,0,...) provided a natural split of the single vector equation into one for a scalar function u (equal the 
former component U 1 ) and a new vector function U (equal the remaining components (U 2 , U 3 , ...)). General 
investigations of systems with one scalar and one vector brought surprisingly a rich set of integrable systems. 
Results in the case A = 2 were analysed by Vladimir Sokolov and Takayuki Tsuchida. What appears to be 
a new solution in this class is the vector generalization 

u t = u xxx + 6uu x -l2(U,U x ), (?) 
U t = -2U XXX - 6uU x . 

of the two-component coupled KdV system proposed by Hirota and Satsuma JH] ■ A Lax pair representation 
has been found by Takayuki Tsuchida. 
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Especially the case of A = 1 with 25 solutions analysed by Takayuki Tsuchida posed a major challenge. 
Not all systems are of interest. Some become triangular by defining (U, U) as a new scalar variable. Others 
are just the result of splitting a single vector equation into scalar + vector equations for the scalar U 1 and 
the vector (U 2 , U 3 , . . . , U ). From the interesting cases just one should be shown here: 

^xxx 

6u 2 u x + u x {U, U) + 2u(U, U x ) + (U, U xx ) + (U x , U x ), 

U t = -2U XXX - 6u xx U - 6u x U x + \2uu x U + 6u 2 U x + (U, U)U X (8) 
-2(U,U X )U. 

Its integrability can be established (Tsuchida) through a change to new variables w and W 

r w = Ux + u 2 + i(u,u), 

\W = U X + 2uU, 

which satisfy the following system: 

j w t = w xxx - 6ww x + 2(W, W x ) 
\w t = -2W XXX + QwW x . 

This system coincides up to a scaling of variables with the coupled Hirota-Satsuma system (JJJ) found in the 
A = 2 investigation above. 

The Miura-type transformation is a generalization of the one for scalar U in fl9\ 120] and the one for 
two-component vector U in 21 . 

In the case when U is a scalar variable, one can set 



(9) 



(10) 



u 



-(q + r), U 



-i(q — r) 



(11) 



2^ ' 2 
and rewrite (jSJ) as a system of symmetrically coupled mKdV equations: 

j Qt = \-\q X x + \r X x + 3(g - r)q x - 2r 3 ] B , 
I r t = [\q X x ~ \r xx - 3(q - r)r x - 2q\. 

This system is identical to (63) in [16 or (3.22) in [17J. It was found in connection with the Kac-Moody Lie 
algebras and written in the Hamiltonian form about twenty years ago (cf. the Cjj case in ^J] or the B2 
case in |2Hj). 

The following table gives an overview of systems for one scalar and one vector found by Crack to have 
higher order symmetries. Details will be discussed in a future contribution. 
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orders of (eq,sym) 


(2,3) 
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(3,5) 


(2,3) 
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# of unknowns (eq,sym,tot) 


5,6,11 
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10,39,49 
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36,164,200 


# of equations 


13 


26 


50 


66 
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386 


149 


313 
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total # of terms 


34 


77 


218 


341 


770 


5000 


1093 


3096 


27695 


av. # of terms/equ. 


2.6 


2.9 


4.3 


5.1 


6.3 


13 


7.3 
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24 


time to formulate 


0.5s 


Is 


5s 


1.8s 


5s 


2m 52s 


8s 


48s 


2h 7m 


time to solve 


0.5s 


0.4s 


6.5s 


29s 


lm 58s 


5h 47m 


29s 


3m 44s 


1 day 


# of solutions 








4 


3 


3 


25 
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A 


scalar: 1/3, vec: 2/3 


scalar: 2/3, vec: 1/3 


orders of (eq,sym) 


(2,3) 


(2,4) 


(3,5) 


(2,3) 


(2,4) 


(3,5) 


# of unknowns (eq,sym,tot) 


10,24,34 


10,54,64 


24,115,139 


13,22,35 


13,66,79 


22,126,148 


# of equations 
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215 


798 


114 


276 


955 


total # of terms 


529 


1462 


12694 


694 


2435 


17385 


av. # of terms/equ. 


5.2 


6.8 


16 


6.1 


8.8 


18 


time to formulate 


3.2s 


13s 


23m 45s 


6.3s 


48s 
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time to solve 
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Table 3. A classification of systems of one scalar + one vector equations. 



5 Classification of Integrable Hamiltonians 

In work done together with Olga V. Efimovskaya quadratic Hamiltonians H have been investigated that 
have the form 

H = (U, AU) + (U, BV) + (V, CV) + (N, U) + (M, V), (12) 

where U = (£^1,^2,^3) and V = (Vx,T^,V3) are three dimensional vectors, A, C are symmetrical matrices, 
B is an arbitrary matrix and N, M are constant vectors. Such Hamiltonians are relevant in the dynamics of 
rigid bodies. 

The equations of motion in the rigid body dynamics are defined by a linear Poisson bracket of the form 

{Y i ,Y j } = c%Y k , i,j,k = l,...,N (13) 
where c^- are some constants. The evolution of dynamic variables Yi, . . . , Y/v is defined by the formula 

where H is the Hamiltonian. 

The skew-symmetricity and the Jacobi identity for the linear Poisson bracket is equivalent to the fact 
that c!a are the structural constants of some Lie algebra. It is known that the Hamiltonian structure of most 
cases of rigid body dynamics can be defined by the linear Poisson brackets 

{U i ,U j }=e ljk U k , {U^Vj} = e ljk V k , {V i ,V j } = 0, (14) 

corresponding to the Lie algebra e(3). For example, two classical problems with a Hamiltonian of form (|T2~j) 
and the Poisson structure (J 14)) are 

1. the Kirchhoff problem (where N = 0, M = in H), 

2. the problem of motion of a massive rigid body around a fixed point. In this case B = C = N = 0. 
The bracket (fT4*j) possesses two Casimir functions: 

Ji = Vl + V 2 2 + V 3 2 , J 2 = UM + U 2 V 2 + U 3 V 3 . (15) 

Therefore to integrate a system on e(3) we need one additional first integral /, functionally independent of 
H, Ji, J 2 . All Hamiltonians ()12p admitting an additional polynomial first integral of first or second degree are 
well known. The main goal of our work is to find all Hamiltonians (112)) that admit an additional polynomial 
integral of degree 3 or 4. 

There exist two kinds of linear transformations of U and V which preserve the Poisson structure (|14j) 
and the form of H. The first kind is defined by 

U = TU, V = TV, (16) 

where T is an arbitrary constant orthogonal matrix. The second kind is defined by 

U = U + SV, (17) 
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where S is an arbitrary antisymmetric matrix. With the transformations (|16|) we reduce the matrix A to 
diagonal form: A = diag(ai,a 2 ,a 3 ). Transformations (|17|) are usually used for a simplification of matrix B. 

The generic case a\ 7^ a 2 7^ 0.3 7^ ai has been investigated in detail and all cases when there exists an 
additional polynomial integral are known. In calculations mentioned below we study the case = a 2 7^ a 3 . 
Using transformation (|17j) we can have 612 = &13 = ^23 — 0. Subtracting multiples of both casimirs enables 
bu — c n — an d gives a Hamiltonian of the form 

H = Uf + Ul + a 3 Ul+ 

2b 21 U 2 V 1 + 2631^3^1 + 2b 32 U 3 V 2 + 2b 22 U 2 V 2 + 2633^3^3+ 

(18) 

2c l2 V l V 2 + 2C 13 U 1 U 3 + C 2 2^2 2 + 2C23V2F3 + C 33 V 3 2 + 

PlUx + p 2 U 2 + p 3 U 3 + qiVi + q 2 V 2 + q 3 V 3 . 

We consider Hamiltonians (fTHj) that have an additional cubic integral. The ansatz for a general first 
integral I of third degree involves 80 terms. Together with the 16 unknown constants in ()18j) the bi-linear 
algebraic system which results from {H, 1} = involves 96 unknowns to be determined. The result of the 
computation is summarized in the following theorem. 

Theorem 1. 

• The Hamiltonian ()18|) admits a polynomial integral of third degree iff it has the form 

H = Uf + Ui + 8l U$ + s 2 V 3 U 3 + s 3 Vi + s A U 3 + s 5 V 3 , (19) 
where Si are arbitrary parameters; 

• The Hamiltonian (fTHj) admits a polynomial integral of third degree on a special level J 2 = of Casimir 
function < li<5j) iff either it has the form hi 9(1 or the form 

H = Uf + Ui + AUi + 4(siVi + s 2 V 2 )U 3 - {s\ + sl)V 3 2 + 

(20) 

s 3 U 3 + S4V1 + S5V2, 

where Sj are arbitrary parameters. 

Hamiltonian (J18)) is a trivial generalization of the Hamiltonian for the Lagrange and Kirchhoff classical 
integrable cases. Actually Hamiltonian (JTHj) admits not only a third degree but also an additional first degree 
integral I = U 3 . If s± — s 2 — then Hamiltonian ()20j) describes the so-called Goryachev-Chaplygin case in 
the problem of motion of a rigid body around a fixed point. The integrability of the general Hamiltonian 
(j!8)l has been recently established in [Tlj . 

6 Non-local 2+1 Dimensional Equations 

The computation decribed in this section solves only a first special case of a wider problem. We still show it 
as it gives an example of how even non-local 2+1 dimensional classification problems can be reduced to the 
solution of bi-linear algebraic systems for which Crack can be useful. 
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The Kadomtsev-Petviashvili equation 



Ufx \^Utl x ~t~ ^xxx)x ~t~ Uyy 

can also be written as 

u t = -(Quu x + u xxx ) + A 2 u x 

with A = D~ 1 D y . A. Mikhailov and R.I. Yamilov observed in ^3] that all known integrable 2+1 - dim 
equations can be written as 

u t = expression in A, D x and u. 

Based on this idea A. Mikhailov, V.V. Sokolov and R. Hernandez Heredero did work at classifying KdV- 
type integro-differential equations of the form 

u t = au xxx + bA(u xxx ) + cA 2 (u xxx ) + ... 

+eA~ s (u xxx ) + terms of lower x-order 

which have symmetries of the same form. 

The special ansatz for equation plus symmetry that has been investigated with Crack is 

u t = u xxx + bA{u xxx ) + cA 2 {u xxx ) 

+0 • A'\u xxx ) + eA" 1 ^^) + lower order (21) 
u T = ■ Uxxx + b'A(u xxx ) + c'A 2 (u xxx ) 

+A'\u xxx ) + e'A~ l (u xxx ) + lower order 

with two coefficients equal one and two zero coefficients due to suitable linear combinations of both equations. 
The lower order terms can be of x-order 2 at most, linear or quadratic in u, and with A restricted as above. 

For this ansatz, which is the simplest in this class, the condition u tT — u Tt = provides already 2865 
bilinear algebraic conditions for the 70 unknowns b,c, . . . and 70 unknowns b', c', . • • ■ 
According to CRACK only 2 solutions exist. One is the Boiti, Leon, Manna, Pempinelli equation 

u t = u xxx + aD x {uA~ l {u)) 

with symmetry 

u T = A 3 (u xxx ) + D x (A(uA(u))) 

The second solution is the same for u — > Au. 

Other known 2-dimensional integrable equations, like KdV in the form u t = Au xxx + 4uAu x + 2u x Au 
which is a symmetry of the usual KdV equation would not have terms u xxx and A 3 (u xxx ) as it was required 
in the special ansatz (J21j) . 

Appendix: A short description of CRACK 

Any identifiers or numbers in curled brackets { } provided at the end or within the following paragraphs refer 
to key-words, file names, module numbers or flags which can be looked up in the Crack manual crack.tex 
(see below under 'Availability') or even be searched in the source code if needed. 
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Philosophy: The program Crack is a computer algebra package written in REDUCE for the solution of 
over-determined systems of algebraic, ordinary or partial differential equations with at most polynomial 
non-linearity. It was originally developed to run automatically and effort was taken for the program to 
decide which computational steps are to be done next with a choice between integrations, separations, 
substitutions and investigation of integrability conditions. It is known from hand computations that 
the right sequence of operations with exactly the right equations at the right time is often crucial to 
avoid an explosion of the length of expressions. This statement keeps its truth for the computerized 
solution of systems of equations as they become more complex. As a consequence more and more 
interactive access has been provided to inspect data, to specify how to proceed with the computation 
and how to control it. This allows the human intervention in critical stages of the computations, {off' 
batch_mode} 

General Structure: A problem consists of a system of equations and a set of inequalities. With each 
equation are associated a short name and numerous data, like size, which functions, derivatives and 
variables occur but also which investigations have already been done with this equation and which not 
in order to avoid unnecessary duplication of work. These data are constantly updated if the equation 
is modified in any way. 

A set of about 30 modules is available to integrate, substitute, decouple, ... equations. A complete list 
can be inspected in interactive mode with the command p2, each operation is listed with the number it 
is called. All modules can be called interactively or automatically. Automatic computation is organized 
by a priority list of modules (each represented by a number) where modules are invoked in the order 
they appear in the priority list, each module trying to find equations in the system it can be applied 
to. If a module is not successful then the next module in the list is tried, if any one is successful then 
execution starts again at the beginning of the priority list. { prog_list_, defaultjprocJ,ist-, fulLprocJist-} 

Because each module has access to all the data, it is enough to call a module by its number. For 
example, the input of the number 2 in interactive mode will start the direct separation module (see 
below) to look for a directly separable equation and will split it. 

Modules: The following modules are represented by numbers in the priority list. Each module can appear 
with modifications under different numbers. For example, integration is available under 7, 24 and 25. 
Here 7 encodes an integrations of short equations = d n f/dx n . 7 has highest priority of the three 
integrations. 24 encodes the integration of an equation that leads to the substitution of a function and 
25 refers to any integration and has lowest priority. 

Integration and Separation: An early feature in the development of the package Crack was the 
ability to integrate exact differential equations and some generalizations of them (see |25j). As 
a consequence of integrations an increasing number of functions of fewer variables is introduced 
which sooner or later produces equations with some independent variables occuring only explicitly 
and not as variables in functions. Such equations are splitted by the integration module. 

Substitutions: Substitutions can have a dramatic effect on the size and complexity of systems. There- 
fore it is possible to have them not only done automatically but also controlled tightly, either by 
specifying exactly which unknown should be substituted where using which equation, or by pick- 
ing a substitution out of a list of substitutions offered by the program {cs} . Substitutions to be 
performed automatically can be controlled with a number of filters, for example, by 

• limiting the size of the equation to be used for substitution, {lengthJimit} 

• limiting the size of equations in which the substitution is to be done, {pdelimit} 
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• allowing only linear equations to be used for substitutions, {linsubst} 

• allowing equations to increase in size only up to some factor in order for a substitution to be 
performed in that equation, {costJimit} 

• allowing a substitution for a function through an expression only if that expression involves 
exclusively functions of fewer variables, {less-vars} 

• allowing substitutions only that do not lead to a case distinction coefficient = or not, 

• specifying whether extra effort should be spent to identify the substitution with the lowest 
bound on growth of the full system. {min_growth} 

Substitution types are represented by different numbers depending on the subset of the above 
filters to be used. If a substitution type is to be done automatically then from all possible 
substitutions passing all filters of this type that substitution is selected that leads to no sub-cases 
(if available) and that uses the shortest equation. 

Factorization: It is very common that big algebraic systems contain equations that can be factorized. 
Factorizing an equation and setting the factors individually to zero simplifies the whole task 
because factors are simpler expressions than the whole equation and set to zero they may lead to 
substitutions and thereby further simplifications. The downside is that if problems with, say 100 
unknowns, need 40 case-distinctions in order to be able to solve automatically for the remaining 
60 unknowns then this would require 2 40 cases to be investigated which is impractical. The 
problem is to find the right balance, between delaying case-distinctions in order not to generate 
too many cases and on the other hand introducing case distinctions as early as necessary in order 
to simplify the system. This simplification may be necessary to solve the system but in any case 
it will speed up its solution (although at the price of having to solve a simplified system at least 
twice, depending on the number of factors). 

For large systems with many factorizable equations the careful selection of the next equation to 
be factorized is important to gain the most from each factorization and to succeed with as few as 
possible factorizations. Criteria which give factors and therefore equations a higher priority are 

• the number of equations in which this factor occurs, 

• if the factor is a single unknown function or constant, then the number of times this unknown 
turns up in the whole system, 

• the total degree of the factor, 

• the number of factors of an equation, 

• and others. 

It also matters in which order the factors are set to zero. For example, the equation = ab 
can be used to split into the 2 cases: 1. a = 0, 2. a ^ 0, b = or to split into the 2 cases 
I. b — 0, 2. b 7^ 0, a — 0. If one of the 2 factors, say b, involves functions which occur only 
linearly then this property is to be preserved and these functions should be substituted as such 
substitutions preserve their linearity. But to have many such substitutions available, it is useful 
to know of many non-linearly occuring functions to be non-zero as they occur as coefficients of the 
linearly occuring functions. In the above situation it is therefore better to do the first splitting 
1. a = 0, 2. a 7^ 0,b = because o ^ will be more useful for substitutions of linear functions 
than b ^ would be. 

An exception of this plausible rule occurs towards the end of all the substitutions of all the 
linearly occuring bi when some b; L are an overall factor to many equations. If one would then 
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set, say, 6 22 = as the second case in a factorization, the first case would generate as subcases 
factorizations of other equations where 622 = would be the second case again and so on. To 
avoid this one should investigate 622 = as the first case in the first factorization. 
The only purpose of that little thought experiment was to show that simple questions, like which 
factored equation should be used first for case-distinctions and in which order to set factors to 
zero can already be difficult to answer in general. 

Elimination (Grobner Basis) Steps: To increase safety and avoid excessive expression swell one 
can apart from the normal call {30} request to do Grobner basis computation steps only if they 
are simplification steps replacing an equation by a shorter equation. {27} 

In a different version only steps are performed in which equations are included which do not 
contain more than 3 unknowns. This helps to focus on steps which are more likely to solve small 
sub-systems with readily available simple results. {57} 

Often the computationally cheapest way to obtain a consistent (involutive) system of equations 
implies to change the ordering during the computation. This is the case when substitutions of 
functions are performed which are not ranked highest in a lexicographical ordering of functions. 
But Crack also offers an interactive way to 

• change the lexicographical ordering of variables, {ov} 

• change the lexicographical ordering of functions, {of} 

• give the differential order of derivatives a higher or lower priority in the total ordering than 
the lexicographical ordering of functions, {og} 

• give either the total differential order of a derivative of a function a higher priority than 
the lexicographical ordering of the derivative of that function or to take the lexicographical 
ordering of derivatives as the only criterium. {of} 

Solution of an under-determined differential equation: When solving an over-determined sys- 
tem of linear differential equations where the general solution involves free functions, then in the 
last computational step often a single equation for more than one function remains to be solved. 
Examples are the computation of symmetries and conservation laws of non-linear differential equa- 
tions which are linearizable. In Crack two procedures are available, one for under-determined 
linear ODEs {22} and one for linear PDEs, {23} both with non-constant coefficients. 

Indirect Separation: Due to integrations new functions of fewer variables are introduced. Substitut- 
ing functions may lead to equations where no function depends on all variables but all variables 
appear as variables to unknown functions, e.g. = f(x)+g(y) although usually much more compli- 
cated with 10 or 20 independent variables and many functions depending on different combinations 
of these variables. Because no variable occurs only explicitly, direct separations mentioned above 
are not possible. Two different algorithms, one for linear indirectly separable equations {10, 26} 
and one for non-linear directly separable equations {48} provide systematic ways of dealing with 
such equations. 

Indirectly separable equations always result when an equation is integrated with respect to dif- 
ferent variables, like = f xy to / = g(x) + h(y) and a function, here f(x,y), is substituted. 

Function and variable transformations: In the interactive mode one can specify a transformation 
of the whole problem in which old functions and variables are expressed in a mix of new functions 
and variables. 

We conclude the listing of modules and continue with other aspects of Crack.. 
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Exploiting Bi-linearity: In bi-linear algebraic problems we have 2 sets of variables a±, ..,a m and b\,..,b n 
such that all equations have the form = YX=i lk a i k ^j k i Ik £ G. Although the problem is linear in the 
cij and linear in the bj it still is a non-linear problem. A guideline which helps keeping the structure of 
the system during computation relatively simple is to preserve the linearity of either the a, or the bj 
as long as possible. In classification problems of integrable systems the ansatz for the symmetry/first 
integral involves usually more terms and therefore more constants (called bj in applications of Crack) 
than the ansatz for the integrable system (with constants a*). A good strategy therefore is to keep the 
system linear in the bj during the computation, i.e. to 

• substitute only a bj in terms of a*, bk, or an a* in terms of an but not an a* in terms of any bk, 

• do elimination steps for any bj or for an a, if the involved equations do not contain any bk, 

The proposed measures are effective not only for algebraic problems but for ODEs/PDEs too (i.e. to 
preserve linearity of a sub-set of functions as long as possible). {flin_} 

Flexible Process Control: Different types of over-determined systems are more or less suited for an au- 
tomatic solution. With the currrent version (2002) it is relatively save to try solving large bi-linear 
algebraic problems automatically. Another well suited area concerns over-determined systems of linear 
PDEs. In contrast, non-linear systems of PDEs most likely require a more tight interactive control. 
Different modes of operation are possible. One can 

• perform one {a} or more computational steps {g} automatically, each step trying modules in the 
order defined by the current priority list {pi} until one module succeeds in its purpose, 

• perform one module a specific number of times or as long as it is successful, {/} 

• set a time limit until which the program should run automatically, {time-limit, limit-time} 

• interrupt an on-going automatic computation and continue the computation interactively, {stop_} 

• arrange that the priority list of modules changes at a certain point in the computation when the 
system of equations has changed its character, 

• induce distinction whether a user-given expression is zero or not. {44} 

Apart from flexible control over what kind of steps to be done, the steps themselves can be controlled 
more or less too, e.g. whether equations are selected by the module or the user. 

Highest priority in the priority list have so-called to-do steps. The list of to-do steps is usually empty 
but can be filled by any successful step if it requires another specific step to follow instantly. For 
example, if a very simple equation = f x is integrated then the substitution of / should follow straight 
away, even if substitutions would have a low priority according to the current priority list. 

Total Data Control: To make wise decisions of how to continue the computation in an interactive session 
one needs tools to inspect large systems of equations. Helpful commands in Crack print 

• equations, inequalities, functions and variables, {e, pi, f} 

• the occurence of all derivatives of selected functions in any equation, {v} 

• a statistics summary of the equations of the system, {s} 

• a matrix display of occurences of unknowns in all equations, {pd} 

• the value of any LISP variable, {pv} 
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• the value of algebraic expressions that can be specified using equation names 
(e.g. coef fn(e_5,df (f ,x,y) ,2)), {pe} 

• not under-determined subsystems, {ss} 

Safety: When working on large problems it may come to a stage where computational steps are necessary, 
like substitution, which are risky in the sense that they may simplify the problem but also complicate 
it by increasing its size. To avoid this risk a few safety features have been implemented. 

• At any time during the computation one can save a backup of the complete current situation in 
a file and also load a backup, {sb, rb} 

• All key strokes are automatically recorded in a list and are available after each interactive step, 
or when the computation has finished. This list can be fed into Crack at the beginning of a 
new computation so that the same operations are performed automatically that were performed 
interactively before. The purpose is to be able to do an interactive exploration first and to repeat 
it afterwards automatically without having to note with pen or pencil all steps that had been 
done, {history-, old-history} 

• During an automatic computation the program might start a computational step which turns 
out to take far too long. It would be better to stop this computation and try something else 
instead. But in computer algebra with lots of global variables involved it is not straight forward 
to stop a computation in the middle of it. If one would use time as a criterion then it could 
happen that time is up during a garbage collection which to stop would be deadly for the session. 
Crack allows to set a limit of garbage collections for any one of those computations that have 
the potential to last forever, like algebraic factorizations of large expressions. With such an 
arrangement an automatic computation can not get stuck anymore due to lengthy factorizations, 
searches for length reductions or elimination steps. {max_gc-elimin, max-gcjac, max-gc-redJen, 
max-gc-short, max-gcss] 

• Due to a recent (April 2002) initiative of Winfried Neun the parallel version of the computer 
algebra system REDUCE has been re-activated and is running on the Beowulf cluster at Brock 
University . This allows conveniently (with a 2- letter command) to duplicate the current status 
of a CRACK computation to another computer, to try out there different operations (e.g. risky 
ones) until a viable way to continue the computation is found without endangering the original 
session, {pp} 

Managing Solutions: Non-linear problems can have many solutions. The number of solutions found by 
Crack can even be higher because to make progress Crack may have factorized an equation and 
considered the two cases a = and a^O whereas solutions in both cases could be merged to only one 
solution without any restriction for a. This merging of solutions can be accomplished with a separate 
program merge_sol() after the computation. 

Another form of post-processing is the production of a web page for each solution, like 
http : //lie .math. brocku. ca/twolf /bl/vll05o35-sl .html . 

If in the solution of over-determined differential equations the program performs integrations of equa- 
tions before the differential Grobner basis was computed then in the final solution there may be 
redundant constants or functions of integration. Redundant constants or functions in a solution are 
not an error but they makes solutions appear unnecessarily complicated. In a postprocessing step these 
functions and constants can be eliminated. {adjust_fnc, drop-Const(), dropredundant()} 
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Parallelization: The availability of a parallel version of Crack was mentioned above allowing to try out 
different ways to continue an ongoing computation. A different possibility to make use of a cluster of 
computers with Crack is, to export automatically the investigation of sub-cases and sub-sub-cases to 
different computers to be solved in parallel. 

It was explained above how factorizations may be necessary to make any progress but also their 
potential of exploding the time requirements. By running the computation on a cluster and being able 
to solve many more cases one can give factorizations a higher priority and capitalize on the benefit of 
factorizations, i.e. the simplification of the problem. 

Relationship to Grobner Basis Algorithms: For systems of equations in which the unknown constants 
or functions turn up only polynomially a well known method is able to check the consistency of the 
system. For algebraic systems this is the Grobner Basis Method and for systems of differential equations 
this is the differential Grobner Basis method. To guarantee the method to terminate a total ordering of 
unknowns and their derivatives has to be introduced. This ordering determines which highest powers 
of unknowns are to be eliminated next or which highest order derivatives have to be eliminated next 
using integrability conditions. Often such eliminations lead to exponential growth of the generated 
equations. In the package Crack such computations are executed with only a low priority. A higher 
priority have operations which reduce the length of equations, irrespective of any orderings. Violating 
any ordering a finite number of times still guarantees a finite algorithm. The potential gain is large as 
described next. 

Length Reduction of Equations: An algorithm designed originally to length-reduce differential equa- 
tions proved to be essential in length reducing systems of bi-linear algebraic equations or homogeneous 
equations which resulted from bi-linear equations during the solution process. 

The aim of the method is to find out whether one equation = E\ can length reduce another one 
= E 2 by replacing E 2 through an appropriate linear combination aE\ — /3E 2 , (3^0. To find a,/3 
one can divide each term of E 2 through each term of E\ and count how often each quotient occurs. If 
a quotient a/ (3 occurs m times then aE\ — (3E 2 will have < n\ + n 2 — 2m terms because 2m terms will 
cancel each other. A length reduction is found if n\ + n 2 — 2m < max(ni, n 2 ). The method becomes 
efficient after a few algorithmic refinements discussed in |24j . Length reduced equations 

• are more likely to length reduce other equations, 

• are much more likely to be factorizable, 

• are more suited for substitutions as the substitution induces less growth of the whole systems and 
introduces fewer new occurences of functions in equations, 

• are more likely to be integrable by being exact or being an ODE if the system consists of differential 
equations, 

• involve on average fewer unknowns and make the whole system more sparse. This sparseness can 
be used to plan better a sequence of eliminations. 

Customization: The addition of new modules to perform new specialized computations is easy. The new 
module only has to accept as input the system of equations, list of inequalities, list of unknowns to 
be computed and provide output in a similar form. The module name has to be added to a list of all 
modules and a one line description has to be added to a list of descriptions. This makes it easy for 
users to add special techniques for the solution of systems with extra structure. A dummy template 
module {58} is already added and has only to be filled with content. 
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Debugging: A feature useful mainly for debugging is that in the middle of an ongoing interactive computa- 
tion the program can be changed by loading a different version of Crack procedures. Thus one could 
advance quickly close to the point in the execution where an error occurs, load a version of the faulty 
procedure that gives extensive output and watch how the fault happens before fixing it. 

The possibility to interrupt REDUCE itself temporarily and to inspect the underlying LISP environ- 
ment {br} or to execute LISP commands and to continue with the Crack session afterwards {pc} 
led to a few improvements and fixes in REDUCE itself. 

Availability: The package Crack including manual can be downloaded for free from 

http://lie.math.brocku.ca/twolf/crack/. It is requested to cite this paper in a publication if 
Crack has been used for any computations that contributed to that publication. 
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